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ABSTRACT 

We use hydrodynamic simulations with realistic models for stellar feedback to study galaxy mergers. 
These high resolution (~ 1 pc) simulations follow the formation and destruction of individual giant molec- 
ular clouds (GMCs) and star clusters. We find that the final starburst is dominated by in situ star formation, 
fueled by gas which flows inwards due to global torques. The resulting high gas density results in rapid 
r j star formation. The gas is self gravitating, and forms massive (< 1O 1O M ) GMCs and subsequently super- 

starclusters (with masses up to 1O 8 M ). However, in contrast to some recent simulations, the bulk of new 
stars which eventually form the central bulge are not born in superclusters which then sink to the center 
Oh of the galaxy. This is because feedback efficiently disperses GMCs after they turn several percent of their 

q mass into stars. In other words, most of the mass that reaches the nucleus does so in the form of gas. 

The Kennicutt-Schmidt law emerges naturally as a consequence of feedback balancing gravitational col- 
^ lapse, independent of the small-scale star formation microphysics. The same mechanisms that drive this 

relation in isolated galaxies, in particular radiation pressure from IR photons, extend, with no fine-tuning, 
over seven decades in star formation rate (SFR) to regulate star formation in the most extreme starburst 
r ^ systems with densities > 1O 4 M0 pc -2 . This feedback also drives super- winds with large mass loss rates; 

^ however, a significant fraction of the wind material falls back onto the disks at later times, leading to 

higher post-starburst SFRs in the presence of stellar feedback. This suggests that strong AGN feedback is 
t-H required to explain the sharp cutoffs in star formation rate that are observed in post-merger galaxies. 

We compare the results to those from simulations with no explicit resolution of GMCs or feedback 
("effective equation of state" [EOS] models). We find that global galaxy properties are similar between 
EOS and resolved-feedback models. The relic structure and mass profile, and the total mass of stars formed 
in the nuclear starburst are quite similar, as is the morphological structure during and after merger (tails, 
(S| bridges, etc.). Disk survival in sufficiently gas-rich mergers is similar in the two cases, and the new models 

t—H follow the same scalings derived for the efficiency of disk re-formation after a merger as derived from 

previous work with the simplified EOS models. While the global galaxy properties are similar between 
EOS and feedback models, sub-galaxy scale properties and the star formation rates can be quite different: 
the more realistic models exhibit significantly higher star formation in tails and bridges (especially in 
shocks), and allow us to resolve the formation of super star-clusters. In the new models, the star formation 
is more strongly time variable and drops more sharply between close passages. The instantaneous burst 
enhancement can be higher or lower, depending on the details of the orbit and initial structural properties 
of the galaxies; first passage bursts are more sensitive to these details than those at final coalescence. 

Key words: galaxies: formation — galaxies: evolution — galaxies: active — star formation: general — 
cosmology: theory 
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1 INTRODUCTION 

A wide range of observed phenomena indicate that gas-rich merg- 
ers are important to galaxy evolution and star formation. In the local 
Universe, the population of star-forming galaxies appears to tran- 
* E-mail :phopkins@ astro.berkeley.edu sition from "quiescent" (non-disturbed) disks (which dominate the 
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total star formation rate/IR luminosity density) at the luminous in- 
frared galaxy (LIRG) threshold lO n L (M* - lO-2OM yr- 1 ) 
to violently disturbed systems at a few times this luminosity. The 
most intense starbursts at z = 0, ultraluminous infrared galaxies 
(ULIRGs; Lir > 10 12 Lq), are invariably associated with mergers 
(e.g.|Joseph & Wright| 1983) | Sanders & Mirabel|1996[|Evans et al.| 
2009 a]), and are fuel ed by compact, central concentrations of gas 
( [Scoville et al.| 1986, Sar gent et al.|1 987) which provide material to 
feed black hole (BH) growth ( [Sanders et al.|1988| ), and to boost the 
concentration and central phase space density of merging spirals to 
match those of ellipticals (Hernqui st et al.|19 93, Robertson et al.| 
2006b). With central densities as large as ~ 1000 times those in 
Milky Way giant molecular clouds (GMCs), these systems provide 
a laboratory for studying star formation under the most extreme 
conditions. 

In addition, various studies have shown that the mass involved 
in these starburst events is critical for explaining the relations be- 
tween spirals, mergers, and ellipticals, and has a dramatic impact 



on the properties of merger remnants (e.g., Lake & Dressier 1986 



Doyon etal.|199 4 Shier & Fischer|l 998 [| James et al.|1999[|G enzel 
et al.|20011|Tacconi et al.|2002||Dasyra et al.|2006||20071|Rot hberg 
& Joseph|2004||2006[|Hopkins et al.|2009a|d) . Even a small mass 
fraction of a few percent formed in these nuclear starbursts can have 
dramatic implications for the mass profile structure ( Mihos & Hern- 
|quist| 1994a} , phase space densities (Hernquist et al.|1993| ), rotation 
and higher-order kine matics (|Cox et al.| 2006b ), kinematically de- 
coupled compo nents ([Hoffman et al.||2009[|20T0]>, s tellar popula- 
tion gradients ( |Soto & Martin||2010[ |Kewley et al.||2010[ |Torrey| 
et al.|2012|), and growth of the central BH (|Di Matteo et al.|2005[ 



Hopkin s et al.|200 6 , Hopkins & Hernquist 2009 ). 

At higher redshifts (z ~ 1 — 3), galaxies with the luminosi- 
ties of LIRGs and ULIRGs may be more "normal" galaxies, in the 
sense that they are relatively undisturbed disks rather than merg- 
ers (Yan et al.|[2007| [Sajina et al.||20Q7| |Dey et al.||2008[ peT] 
|bourne et al. 2008 ; Dasy ra et al.||2008| L However, at those same 
redshifts, yet more luminous systems appear, including large popu- 
lations of Hyper-LIRGs (Hy LIRGs; Lir > 1O 13 L ) and bright sub- 
millimeter galaxies ( e.g.|Chapman et al.|2005[| Younger et aTj 2007 
|2009b[ |Casey et al.||2009| ). These hyperluminous objects exhibit 
many of the traits associated with merger-driven starbursts, includ- 
ing morphological disturbances, and may be linked to the emer- 
gence of massive, quenched (non star-forming), compact ellipticals 
at times as early as z ~ 2 — 4 ( Papovich et al. 2005 ; Younger et al. 
2008 b[|Tacconi et al.|2006l |Schinnerer et al.|2008[|Chapman et al." 



2009, Tacconi et al.|2008 ). Reproducing their abundance and lumi- 



nosities remains a challenge for current models of galaxy formation 
([Baugh et al.|2005[l Swinban k et al.|2008l [Narayanan et al.|2010[ 
Youn ger et al'|2009a[|Hayward et al.poTo] |2011a| ). 

Modeling the consequences of these mergers for star for- 
mation and the nuclear structure of galaxies requires following 
the highly non-linear, resonant, and chaotic interplay between gas 
(with shocks, cooling, star formation, and feedback), stars, and dark 
matter over a very large dynamic range. High-resolution numerical 
hydrodynamic simulations are the method of choice. However, un- 
til recently, computational limitations have made it impossible to 
resolve the ^pc (and < lOOOyr) scales of structure in the ISM rele- 
vant for star formation and stellar feedback while following the four 
or five orders -of-magnitude larger global evolution of a galaxy or 
merger. Perforce, most previous models have adopted a variety of 
sub-grid approaches to describe "effective" ISM properties below 
the resolution scales of the simulations. 

All numerical galaxy simulations require some sort of "sub- 



grid" model; for example, even with box sizes of order one par- 
sec, current star formation simulations treat radiative feedback in 
a very approximate manner, while it is clear that such feedback 
is crucial to determining, e.g., the inital mass function. Unfortu- 
nately, without a physically well founded sub-grid model model, 
the predictive power of numerical simulations is limited. If the av- 
erage star formation properties are put in by hand (and they are in 
most galaxy scale simultions to date) the global star formation rate 
clearly cannot be predicted. Similarly, if the effects of stellar feed- 
back are put in by hand, e.g., by invoking an "effective pressure" 
or turbulent dispersion, in order to suppress runaway cooling and 
clumping, the detailed physics and implications for star formation 
cannot be determined. Moreover, it is by no means clear if com- 
monly used prescriptions capture the key physics, nor whether they 
can be extrapolated from "quiescent" systems (which are typically 
used to calibrate the models) to the very different ISM conditions 
in mergers. Typically subgrid models also have several adjustable 
parameters, which further limit their predictive power; it is often 
unclear whether observational constraints imply differences in the 
merger dynamics or more subtle tweaks in the sub-grid ISM as- 
sumptions. This has occasionally led to contradictory conclusions 
in the literature. 

Numerical simulations of isolated galaxies and galaxy merg- 
ers can now achieve the dynamic range required to resolve the for- 
mation of GMCs and ISM structure, ~ 1 — 10 pc (se e e.g. [S aitoh 
et al.|2008[|Tasker & Tan|2009l|Bournaud et al.|2010a||Dobbs et all 



2011| ). However, improving resolution alone - beyond the implicit 
"averaging" scales for the ISM model - has no clear meaning if 
the physics that govern star formation and ISM structure on these 
scales are not included. Generally, models have not attempted to 
include these physics or have included a very limited sub- set of 
the relevant processes. In fact, a large number of feedback mech- 
anisms may drive turbulence in the ISM and help disrupt GMCs, 
including: photo-ionization, stellar winds, radiation pressure from 
UV and IR photons, proto- stellar jets, cosmic rays, supernovae, and 
gravitational cascades from large scales (e.g. |Mac Low" & Klessen 
|2004[ and references therein). 

In |Hopkins et"aL] ( |2011b| ) (Paper I) and |Hopldns"eTaTl ( [20T2l ) 
(Paper II) we developed a new set of numerical models to incorpo- 
rate feedback on small scales in GMCs and star-forming regions, 
into simulations with pc-scale resolution^] These simulations in- 
clude the momentum imparted locally (on sub-GMC scales) from 
stellar radiation pressure, radiation pressure on larger scales via 
the light that escapes star-forming regions, HII photoionization 
heating, as well as the heating, momentum deposition, and mass 
loss by SNe (Type-I and Type-II) and stellar winds (from O and 
AGB stars). The feedback is tied to the young stars, with the en- 
ergetics and time-dependence taken directly from stellar evolution 
models. Our method also includes realistic cooling to temperatures 
< 100K, and a treatment of the molecular/atomic transition in gas 
and its effect on star formation. 

We showed in Papers I & II that in isolated disk galaxies these 
feedback mechanisms produce a quasi- steady ISM in which giant 
molecular clouds form and disperse rapidly, after turning just a few 
percent of their mass into stars. This leads to an ISM with phase 
structure, turbulent velocity dispersions, scale heights, and GMC 
properties (mass functions, sizes, scaling laws) in reasonable agree- 
ment with observations. In Hopkins et al. (2011c ) (Paper III), we 



1 Movies of these simulations are available at |https : //www, cf a . | 
|harvard . edu/ ~phopkins/Site/Movies_sbw_mgr . html 
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showed that these same models of stellar feedback predict the elu- 
sive winds invoked in almost all galaxy formation models; the com- 
bination of multiple feedback mechanisms is critical to give rise to 
massive, multi-phase winds having a broad distribution of veloci- 
ties, with material both stirred in local fountains and unbound from 
the disk. 

Papers I & II showed that the global star formation rate found 
in the simulations did not depend on the sub-grid model for star 
formation, over a broad range of sub-grid model parameters and 
even model types. The conclusion drawn was that the rate of star 
formation in the simulations was controlled by the amount of feed- 
back required to maintain the gas disk in hydrostatic equilibrium, 
and with a Toomre Q parameter of order unity. If this is true in real 
galaxies, it will simplify the task of understanding galaxy formation 
greatly. 

In this paper, we extend these models to idealized major merg- 
ers between galaxies. In this first exploration, we compare some of 
the global properties of the remnants - particularly those tied to 
star formation and ISM physics - to the results of simulations us- 
ing more simplified sub-grid ISM/feedback models. Specifically, 
we investigate the consequences of a more realistic treatment of 
the ISM for the star formation histories, spatial distribution of stel- 
lar mass and starburst stars, survival and re-formation of gas disks, 
and origins of the global Kennicutt- Schmidt law. All these proper- 
ties can be compared to the results from previous sub-grid models 
in a straightforward manner. We also test the proposition that the 
star formation rate is controlled by feedback, in mergers as well 
as in quiescent galaxies, since we do not adjust our star formation 
prescription. 

In companion papers, we examine in more detail some of 
the phenomena that could not be predicted with previous models: 
the phase structure of the ISM and properties of starburst "super- 
winds" driven by stellar feedback, as well as the physics of star 
cluster formation in tidally shocked and starburst regions. 



t = 0.5 
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Figure 1. Morphology of the gas in a standard simulation (high-resolution 
with all explicit feedback mechanisms included) of a merger of the HiZ 
disk model: a massive, z ~ 2 — 4 starburst disk merger. The time is near 
apocenter after first passage. This color projection emphasizes the cold, 
star-forming gas. Brightness encodes projected gas density (light-to-dark 
with increasing density; logarithmically scaled with a « 4dex stretch); color 
encodes gas temperature with the blue/violet material being T < 1000K 
molecular and atomic gas, pink/red ~ 10 4 — 10 5 K warm ionized gas, and 
yellow/white > 10 6 K hot gas. Gravitational collapse forms giant molecular 
clouds and proto-star cluster complexes throughout the gas. The outflows 
present include a component in dense clumps. Images of the other simula- 
tions at various times are in Appendix |A| 



2 METHODS 

The simulation techniques and galaxy models used here are de- 
scribed in detail in Paper I (Sec. 2 & Tables 1-3) and Paper II 
(Sec. 2). We briefly summarize the most important properties here. 
The simulations were performed with the parallel TreeSPH code 
GADGET-3 ([Sp ringel 2005 ), in its fully conservative formulation 
(Springel & Hernquist 2002 ). They include stars, dark matter, and 
gas, with cooling, shocks, star formation, and stellar feedback. 



2.1 Disk Models 

We implement the ISM prescription in four distinct initial disk 
models spanning a range of galaxy types (summarized in Table [TJ. 
Each has a bulge, stellar and gaseous disk, halo, and central BH 
(although to isolate the role of stellar feedback, models for BH 
growth and feedback are disabled). At our standard resolution, each 
model has ^0.3 — 1 x 10 8 total particles, giving particle masses of 
500 — 1 000 Mq and 1 — 5 pc smoothing lengths. A couple of ultra- 
high resolution runs (of isolated versions of these disks) for con- 



vergence tests employ w 10 9 particles with sub-pc resolutionj^The 
disk models include: 

(1) SMC: an SMC-like dwarf, with baryonic mass Mbar = 
8.9 x 10 8 M^ and halo ma ss M ha io = 2 x 10 10 Mq (concentration 
c = 15), a Hernquist] ( |l990| > profile bulge with a mass = 10 7 Mq , 
and exponential stellar {ma = 1.3 x 10 8 Mq) and gas disks (m g — 
7.5 x 10 8 Mq) with scale-lengths hd = 0.7 and h g — 2.1 kpc, re- 
spectively. The initial stellar scale-height is zo = 140 pc and the stel- 
lar disk is initialized such that the Toomre Q—\ everywhere. The 
gas and stars are initialized with uniform metallicity Z = O.IZq. 

(2) MW: a MW-like galaxy, with halo and bary- 
onic properties of (Mhaio, c) — (1.6 x 10 12 Mq, 12) and 
(M har ,m b ,m d ,m g ) (7.1,1.5,4.7,0.9) x 10 10 Mq, Z Z , 
and scale-lengths (hd, h g , zo) = (3.0, 6.0, 0.3) kpc. 

(3) Sbc: a LIRG-like galaxy (i.e. a more gas-rich spiral than is 
characteristic of those observed at low redshifts) with (Mhaio, c) — 
(1.5 x lO n M , 11), (M har ,m b ,m d ,m g ) (10.5, 1.0,4.0,5.5) x 
10 9 M , Z = 0.3Zq, and (h d , h g , zo) = (1.3, 2.6, 0.13) kpc. 

(4) HiZ: a high-redshift massive starburst disk, chosen to 



2 These tests are described in Paper I and Paper II, and used to check con- 
vergence in small-scale ISM properties. We have also run every simulation 
described in this paper with 10 times fewer particles (2 times larger soften- 
ing); although the small-scale properties such as the structure of individual 
star clusters differ, all properties plotted in this paper (more global values 
in general) are found to be indistinguishable. 
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Table 1. Galaxy Models 
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Parameters describing our (isolated) galaxy models, used as the initial conditions for all of the mergers: 

(1) Model name: shorthand for models of an isolated SMC-mass dwarf (SMC); local gas-rich galaxy (Sbc); MW-analogue (MW); and 
high-redshift massive starburst (HiZ). (2) e g : gravitational force softening. Higher-resolution tests of the isolated galaxies use half this value 
(Paper I). (3) Mh a i : halo mass. (4) c: halo concentration. Values lie on the halo mass-concentration relation at each redshift (z = for SMC, 
Sbc, and MW; z = 2 for HiZ). (5) V max : halo maximum circular velocity. (6) M^ ary : total baryonic mass. (7) Aft,: bulge mass. (8) a: 
Hernquist (1990) profile scale-length for bulge. (9) M^: stellar disk mass. (10) r^: stellar disk scale length. (11) h: stellar disk scale-height. 
(12) M g \ gas disk mass. (13) r g : gas disk scale length (gas scale-height determined so that Q = 1). (14) / gas : average gas fraction of the disk, 
inside of the stellar R e (Af g [< R e ]/(M g [< R e ] + Af d [< R e ])). The total gas fraction, including the extended disk, is ~ 50% larger. (15) Z: 
initial metallicity (in solar units) of the gas and stars. 




Figure 2. The stars at the same time as Fig. ^ The image is a mock ugr 
(SDSS-band) composite, with the spectrum of all stars calculated from 
their known age and metallicity, and dust extinction/reddening accounted 
for from the line-of-sight dust mass. The brightness follows a logarithmic 
scale with a stretch of « 2 dex. Young star clusters are visible throughout the 
system as bright while pixels. The nuclei contain most of the star formation, 
evident in their saturated brightness. Fine structure in the dust gives rise to 
complicated filaments, dust lanes, and patchy obscuration of star-forming 
regions. A few super-starclusters are apparent as the brightest young stel- 
lar concentrations. Images of the other simulations at various times are in 
Appendix [A] 

match the properties of the observed non-merging but rapidly star- 
forming SMG population, with (M hak) , c) = (1.4 x 1O 12 M , 3.5) 
and a virial radius appropriately rescaled for a halo at z = 2 rather 
than z = 0, (M bar , m b , m d , m g ) (10.7, 0.7, 3, 7) x 1O 1O M , Z 
0.5Zq, and (hd, h g , zo) = (1.6, 3.2, 0.32) kpc. 

2.2 Merger Parameters 

The purpose of this paper is to contrast models with and without 
explicit stellar feedback, not to present a systematic study of all 



merger parameters - extensive studies of this nature can be found in 
the literature (Barnes & Hernquist 1992, Barnes 2002; Naab et al. 
T999l|Cox et al.||2006b| |2008[ |Robertson et al'||2006b| |Younger 
et al. 2008al|di Matteo et al.|20071|Burkert et al |20081 [Hof fman 
etaL 2009, Jesseit et al. 2009). We therefore focus on a small but 
representative subset of possible mergers. 

We consider equal-mass mergers (merging identical copies of 
galaxies (l)-(4)). The pairs are placed on parabolic orbits (moti- 
vated by cosmological simulations; see Benson 2005 , Khochfar & 
Burkert 2006 ) with the spin axis of each disk specified by and 
4> in spherical coordinates. We choose two representative orienta- 
tions from among those which have been studied in a number of 
other works (|Cox et al.|2006bl|Barnes[T988l|Hernquist|1992||1993l 
|Naab & Burkert|2003| ). The first (orbit e in |Cox et al.|2006b> is near- 
prograde with (Q u <j> u 2 , 4>i ) = (30, 60, -30, 45); the second (or- 
bit f) is near-retrograde (or polar-retrograde) with (6\ , 0i , #2, 4>i) = 
(60,60,150,0). This spans both a strong resonant interaction 
(prograde) and weak, out-of-resonance (retrograde) interaction. 
We choose these rather than perfectly prograde (#1, 0i, 62, 4>i) — 
(0, 0, 0, 0) and retrograde (6>i, </> u 2 , <jn) = (180, 0, 180, 0) orbits 
because there are known pathological behaviors associated with the 
narrow phase space of perfectly resonant orbits. 

We have compared a subset of low-resolution simulations that 
vary orbital inclinations and pericentric passage distances; the vari- 
ation in remnant properties is consistent with previous studies (e.g. 
|Cox et al.|2006b| ). But this variation is largely a dependence on the 
strength of gravitational torques and total angular momentum; the 
differences between models with/without feedback do not dramat- 
ically depend on these properties. In any case our e and f orbits 
reasonably bracket the range from most to least violent orbits. 

2.3 Explicit Feedback Models 

The most important physics in these simulations is the model of 
stellar feedback. We include feedback from a variety of mecha- 
nisms, each of which we briefly describe below. More details about 
our implementations of this physics are given in Paper I and Paper 
II. We use a|Kroupa (2002} initial mass function (IMF) through- 
out and use STARBURST99 ( |Leitherer et al.|1999| > to calculate the 
stellar luminosity, mass return from stellar winds, supernova rate, 
etc. as a function of the age and metallicity of each star particle. 

(1) Local Momentum Deposition from Radiation Pressure, 
Supernovae, & Stellar Winds: In Paper I, we present the radiation 
pressure aspect of this model for feedback from young star clusters 
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Figure 3. Left: Gas, as Fig.[T] at a slightly later time, with a color projection here chosen to emphasize the ionized and hot gas. Brightness now increases with 
surface density; color encodes temperature in a similar manner (blue/pink/yellow representing cold/warm/hot phases). Violent outflows are present, emerging 
both from the complexes and the system as a whole, driven by the massive starburst. The volume-filling "hot" component is now visible, with multiple 
"bubbles" driven by separate local events. Right: Same, but with the (projected in-plane) velocity vectors plotted. The vectors interpolate the gas velocities 
evenly over the image; their length is linearly proportional to the magnitude of the local velocity with the longest plotted corresponding to w 500km s _1 . 






t = 2.5 






♦ 



Figure 4. Morphology of the f (retrograde) merger of the Sbc galaxy (a dwarf starburst); each shows the optical as Fig.|2]at different times during the merger. 



in detail. At each timestep, gas particles identify the nearest density 
peak representing the center of the nearest star-forming "clump" or 
GMC-analog. We calculate the total luminosity of the star particles 
inside the sphere defined by the distance from the center of this 
star-forming region to the gas particle of interest; the incident flux 



on the gas is then determined assuming that the local star forming 
region is optically thick to the UV radiation. 

The rate of momentum deposition from radiation pressure is 
then P r ad ~ (1 + riR)Lincident/c where the term 1 + tir accounts for 
the fact that most of the initial optical/UV radiation is absorbed and 
re-radiated in the IR; tjr = H gas k\ R is the optical depth in the IR, 
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which allows for the fact that the momentum is boosted by multiple 
scatterings in optically thick regions. Here S gas is calculated self- 
consistently as the average surface density of the identified clump, 
with kir « 5 (Z/Zq) g _1 cm 2 approximately constant over the rel- 
evant physical range of dust temperatures. The imparted accelera- 
tion is directed along the flux vector. In Paper I we discuss numer- 
ous technical aspects of this implementation - such as the effects of 
resolution, photon leakage, and how the momentum is discretized 
- and show that essentially all our conclusions are robust to uncer- 
tainties in these choices. 

The direct momentum of SNe ejecta and stellar winds PsNe 
and P w are similarly tabulated from STARBURST99 and injected 
as an appropriate function of age and metallicity to the gas within 
a smoothing length of each star. This source of turbulent energy is 
almost always smaller than that due to radiation pressure discussed 
above. In some cases, however, in particular in dwarf galaxies, the 
work done by bubbles of gas shock-heated by supernovae and/or 
stellar winds is dynamically important; this is discussed below. 

(2) Supernova and Stellar Wind Shock-Heating: The gas 
shocked by supernovae and stellar winds can be heated to high tem- 
peratures, generating bubbles and filaments of hot gas. We tabulate 
the Type-I and Type-II SNe rates from|Mannucci et al. (2006| and 
STARBURST99, respectively, as a function of age and metallicity 
for all star particles, and stochastically determine at each timestep 
if a SN occurs. The SNe are resolved discretely in time (as opposed 
to continuous energy injection). For each SN, the appropriate ther- 
mal energy is injected into the gas within a smoothing length of 
the star particle. Similarly, stellar winds are assumed to shock lo- 
cally and so we inject the appropriate tabulated mechanical power 
L(f , Z) as a continuous function of age and metallicity into the gas 
within a smoothing length of the star particles. The specific energy 
of these stellar winds is large for young stellar populations in which 
fast winds from massive stars dominate, but declines rapidly at later 
times when slower, AGB winds dominate. 

(3) Gas Recycling: Gas mass is returned continuously to the 
ISM from stellar evolution, at a rate tabulated from SNe and stellar 
mass loss in STARBURST99. The integrated mass fraction recy- 
cled is ~ 0.3. 

(4) Photo-Heating of HII Regions: For each star particle, we 
tabulate the rate of production of ionizing photons; starting from the 
nearest gas particle and moving radially outwards, we then ionize 
each non-ionized gas particle (using the gas and stellar properties 
to determine the photon production rate needed to maintain the par- 
ticle as fully ionized). Gas which is ionized is immediately heated 
to rsj 10 4 K, unless is already above this temperature; moreover, the 
gas is not allowed to cool below 10 4 K until it is no longer in an 
HII region. This method allows for overlapping, non-spherical HII 
regions that can extend to radii ^kpc. 

(5) Long-Range Radiation Pressure: Radiation pressure 
from photons absorbed in the immediate vicinity of stars is cap- 
tured in mechanism (1). However, photons that escape these re- 
gions can still be absorbed at larger radii. For each star particle, we 
construct the intrinsic SED (L u ) as a function of age and metallic- 
ity; we then use the local density and density gradients to estimate 
the integrated column density and attenuation of the SED using 

p(h S m\ + | Vlnp| where h sm \ is the smoothing 
length and k v is the frequency-dependent opacity (assuming dust 
opacities that scale with metallicity, as in (1) above). The resulting 
"escaped" SED gives a frequency-dependent flux that is prop- 
agated to large distances. We construct a force tree for this long 
range force in an identical fashion to the gravity tree, since after 
attenuating the flux near the star particle, the stellar flux is assumed 



to decrease oc r -2 . Each gas particle then sees an incident net flux 
vector ¥ l u , integrated over all stars in the galaxy. 

Extensive numerical tests of the feedback models are pre- 
sented in Paper I, Paper II, and Paper III. 

In these models, gas follows an atomic cooling curve with ad- 
ditional fine- structure cooling to < 100K, with no cooling floor 
imposed. Star formation follows the model in Hop kins et al.| 
(2011a). Regions which are locally self-gravitating on the small- 
est resolved scale, i.e. which have a = 5v 2 Sr/Gm ga s(< Sr) —> 
(|V-v| 2 + |Vxv| 2 )/(2Gp) < 1 should collapse to much higher 
densities in a free-fall time absent feedback; we therefore assign 
them an instantaneous SFR of p* = p/fe(p). We stress that be- 
cause feedback is present and can regulate against further collapse 
once stars form, the average efficiency in the dense gas is much 
lower (typically ~ 1%). We further follow Krumholz & Gnedin 
P011| ) and calculate the molecular fraction within the dense gas 
as a function of the local column density and metallicity, and allow 
star formation only from that gas; however as shown in Paper II and 
Hopkins et al. ( 2011a ) this has essentially no effect on our results. 

As noted in the introduction, we showed in Paper I and Paper 
II that the galaxy structure and SFR are basically independent of 
the small-scale SF law, density threshold (provided it is high), and 
treatment of molecular chemistry, because the structure and SFR 
are feedback-regulated. As a further test of this result, we have 
re-run several of the models described here with the more sim- 
plified (and frequently used) star formation prescription in which 
star formation is restricted to occur only in gas above a density 
n > 1000cm -3 , with a rate p* = 0.015p/fe corresponding to ob- 
served average efficiencies "enforced." This has almost no effect 
on global properties or the galaxy- average SFRs, but it does tend to 
"smear out" the SFR in dense regions such as the merging nuclei, 
where all the gas is well above this threshold. 

2.4 Models with an "Effective" Equation of State 

In previous work on galaxy mergers, it was not possible to ex- 
plicitly resolve the processes described above. Instead, simulations 
used a variety of sub-grid approaches to model some "effective" 
ISM properties on the resolution scales of the simulations (often a 
few hundred pc). 

We want to assess the accuracy of such subgrid models. We 
therefore compare otherwise identical disk merger models using 
the effective equation of state (EOS) approach in Spring eT& Hern-| 
|quist| ( |2003| ). Since high-density GMCs are not resolved, the SF 
threshold is lower, chosen to represent densities where multi-phase 
structure should appear (here we adopt n > 1 cm -3 ); the efficiency 
must be correspondingly tuned to match the global Kennicutt re- 
lation (here p* = 0.03 p/ts). The model attempts to account for 
the complex sub-grid feedback physics described above by simply 
assigning the star-forming gas an effective sound speed c e & repre- 
senting e.g. a sub-grid turbulent ISM pressure. 

In the model of Springel & Hernquist (2003 ), the effective 
pressure is motivated by the ISM model in McKee & Ostriker| 
( |1977| ). Because the assumed sub-grid ISM properties are uncer- 
tain, the model allows for interpolation in the EOS with the pa- 
rameter ^eos (see Springel et al. 2005 ; Springel & Hernquist 2005 ; 
Robertson et al.|2006a|c| ). The standard "maximal feedback" case 
(#eos = 1) assumes 100% of the energy from supernovae couples 
to an inter-cloud diffuse phase with long cooling times; this leads 
to a fairly "stiff" equation of state, c e ff ^ 50 — 200km s -1 . The op- 
posite "weak feedback" case (q e0 s = 0) has c e ff °° lOkrns -1 , moti- 
vated by the minimum dispersion observed in gas in most galax- 
ies. This is quite similar to the effective EOS used in |Teyssier| 
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3 MORPHOLOGIES 
3.1 Merger Morphologies 
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Figure 5. Merger remnant stellar mass profiles, with explicit (thick lines) 
and implicit/sub-resolution (thin lines) stellar feedback prescriptions. Top: 
High-redshift gas -rich starburst merger. Middle: MW-like galaxy merger. 
Bottom: SMC-mass dwarf merger. For each simulation, we plot the to- 
tal stellar surface density (black solid), and the separate dissipation- 
less/violently relaxed component (stars formed before the merger simu- 
lation; red dotted) and dissipation/starburst component (stars formed in 
the merger-induced starburst; blue dashed). The profiles are averaged over 
^100 viewing angles. In all cases, the profiles agree very well indepen- 
dent of the details of the feedback prescription. Despite the inner structure, 
dynamics, and winds associated with the starburst being modified consid- 
erably, the "total starburst mass" and characteristic radii where the post- 
starburst light dominates the profile appear robust. This is expected if these 
properties are set by large-scale gravitational torques. 



Figures |1|2| show the gas and stars in some representative stages 
of the merger simulations, for our standard case ("explicit" feed- 
back model) where all feedback mechanisms are present. Figure]?] 
and Figures |A1|A4| (in Appendix [A) extend this to show the gas 
and stars in a number of merger stages, from different viewing an- 
gles, for each merger. The gas maps show the projected gas den- 
sity (intensity) and temperature (color, with blue representing cold 
molecular gas at T < 1000K, pink representing the warm ionized 
gas at ~ 10 4 — 10 5 K, and yellow representing the hot, X-ray emit- 
ting gas at > lO 6 K)0The stellar maps show a mock three-color 
observed image, specifically a u/g/r composite. The stellar lumi- 
nosity in each band is calculated from each star particle according 
to the STARBURST99 model given its age, mass, and metallicity 
(and smoothed over the appropriate kernel). We then attenuate the 
stars following the method of Hopkins et al. (2005): we calculate 
the total dust column (from the simulated gas) along the line-of- 
sight to each star particle for the chosen viewing angle (assuming 
a constant dust-to-metals ratio, i.e. dust-to-gas equal to the MW 
value times Z/Zq), and apply a MW-like extinction and reddening 
curve (as tabulated in |Pei|l9 92)| 4 l 

As seen in Paper II for isolated galaxies, gravitational collapse 
forms GMC complexes in the gas, which are individually dispersed 
rapidly by feedback once they turn a few percent of their mass 
into stars, but are replaced by continuously forming new clumps. 
The ISM as a whole is a a highly supersonically turbulent multi- 
phase medium, with volume-filling hot gas (heated by SNe and O- 
star winds) making up a few percent of the mass, a "warm" phase 
of gas maintained by HII photoheating at ~ 10 4 K, which consti- 
tutes a large filling factor and mass fraction at large radii from the 
galaxy nuclei, and GMC complexes in dense molecular filaments 
(and volume-filling molecular gas near the galaxy centers), with 
most of the gas mass. 

Violent outflows are clearly evident, both emerging from in- 
dividual GMCs and from the merging systems as a whole. The 
composition of these outflows is discussed in detail in Paper III. 
Briefly, the volume filling hot gas is primarily venting SNe and O- 
star wind heated gas; there is some contribution from gravitational 
shocks, but it is much smaller (see |Cox et al.|2006a| and discussion 
below). Occasionally the early stages of these outflows are evident 
as "bubbles" breaking out of the warm phase gas. The warm/cold 
clumps embedded in the outflow predominantly come from mate- 
rial directly accelerated by radiation pressure (which is much more 
efficiently at launching cold gas out of the disk, as opposed to the 
entrainment of that cold material by hot, tenuous outflows). 

In the stellar images, young stars clusters are visible through- 
out the galaxies. Most of the stars form in resolved clusters; in a 
companion paper we discuss the nature and properties of these clus- 
ters in detail. Most of them, however, are open clusters that corre- 
spond to GMCs which have converted ~ 5 — 10% of their mass 
into stars (see Paper II), and so are unbound once the much larger 
(marginally bound) GMC gas mass is dispersed by feedback. The 
bright nuclei contain most of the star formation, as discussed below. 



et al. ( 2010) . Intermediate v alues interpolate as cL = q&o scl^[q 



1] + (1 — g e0 s) c eS [q — 0]. In Hopkins & Quataert (20101), various 



choices are compared to observations of star-forming systems as 
a function of density; they suggest that q CO s ^ 0.1 — 0.3 gives a 
reasonable representation of the observations. We therefore adopt 
g e0 s = 0.25 as our standard choice; but we have also re-run compar- 
ison models with q eos = 1 and our conclusions are nearly identical. 
© 0000 RAS, MNRAS 000, 000-000 



3 The projected temperatures are logarithmically-averaged and surface- 
density weighted, so reflect the temperature of most of the line-of-sight gas 
mass, rather than the temperature that contains most of the thermal energy. 

4 Accounting for dust scattering changes the detailed spectrum, but the col- 



i and morphologies in the images are not significantly altered ( Jonsson 
|et al.|2006||Wuyts et al.|2009||20i0} . 
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The massive gas inflows give rise to large dust column densities 
along many sightlines; face-on, we see heavily extincted dust lanes 
with complex structure (filaments, feathering, etc.). The patchy ob- 
scuration is typical of local merging systems. A few super-star clus- 
ters are visible as the brightest young stellar "knots," often outside 
the nuclei in tidally shocked arms. 

3.2 Merger Remnant Morphologies 

As expected for the remnants of very major (equal-mass) mergers, 
the relic stars typically form a bulge-dominated galaxy once the 
merger has completed. In Figure [5] we show this quantitatively by 
comparing the surface stellar mass density profiles of the remnant 
stars. We plot the median profile averaged over 1000 lines-of- sight 
uniformly sampling the unit sphere, but the sightline-to-sightline 
differences are small. Despite being initially disk-dominated, the 
HiZ and MW models form clear r 1/4 -law bulges; the SMC and Sbc 
cases may still have significant disks (discussed below), evident in 
the curvature (lower Sersic n s ) here. 

The final mass profiles of models with explicit feedback or 
with an effective EOS are quite similar. This is, of course, ex- 
pected where the profile is dominated by violent relaxation of pre- 
existing stars. However, in previous EOS models, it has generally 
been found that the central ^kpc is dominated by stars formed in 
situ in the nuclear starburst (from gas driven to the center by strong 
torques; see |Mihos & Hernquist|1994a[ [Hopkins et al.|2 008b). We 
therefore examine the massprofile of two sub-components: the 
stars formed in the starburstjjand the violently relaxed "envelope" 
of stars formed before this starburst/coalescence (all stars formed 
previous to this cut). In all cases, the explicit feedback models are 
nearly identical to the "effective" equation of state models. In par- 
ticular, the starburst component, not just the dissipationless enve- 
lope, appears to have similar mass, shape, and radius. 

There are some systematic second-order differences. The full 
feedback models tend to produce a somewhat more massive ^kpc- 
scale core (at the tens of percent level), because dissipation is more 
efficient. There are subtle differences at large radii (~ lOkpc) ow- 
ing to extended star formation (discussed below), although the total 
stellar mass in young stars involved at these radii is very small. 

3.3 Disk Survival 

A controversial subject in recent years has been the question of 
disk "survival" or "re-formation" after mergers. Simulations with 
EOS models predicted that sufficiently gas-rich mergers should 
leave large relic disks and even disk-domin ated galaxies ([Spr ingel 



SMC 



& Hernq uist|2005[|Robertson et al.|2006al|Hopkins et al.]20 09b) 
Hop kins et al.| ( 2009b) showed that in those models, this is gov- 
erned by global gravitational torques being inefficient at removing 
angular momentum from gas at large radii in disks, and so is inde- 
pendent of the details of the EOS implementation or introduction 
of an ad-hoc stellar wind model. However, other simulations, with 
either no hydrodynamics ("sticky particle" codes) or weak stellar 
feedback but efficient cooling, predict much less efficient disk sur- 
vival (Bournaud et al. 2010b ). We therefore examine this question 
for the first time in simulations with realistic stellar feedback mod- 
els. 

Figure [6] shows the morphology of one merger remnant (the 
SMC-SMC retrograde f merger) which clearly re-forms a promi- 
nent, rotationally supported disk. We quantify the disk mass and 

5 For simplicity, we define this as all stars formed within ±150Myr of the 
peak in the SFR near coalescence, but small changes in the definition have 
no effect on our conclusions. 
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Figure 6. Disk survival in mergers. Top: Gas image of the SMC f merger 
remnant (as Fig.[TJ, with our full feedback model, after final coalescence. 
The gas has clearly re-formed a large, rotating disk (face-on here), very sim- 
ilar to the isolated progenitor(s). Middle: Distribution of specific angular 
momentum j of the baryons within 5 kpc of the center of the same model, 
after the merger; values are compared to the value j circ of a pure circular 
orbit at the same radius. The distribution is clearly bimodal, with a peak 
containing « 30% of the mass near j/ j c [ rc = (bulge), and « 70% in a 
peak near j/ = 1 (disk). The disk lies on the same Tully-Fisher relation 
as its progenitors. Bottom: Mass fraction of the disks which is consumed in 
the final coalescence starburst (/burst), as a function of the gas mass fraction 
just before the burst begins (/ gas ; here 200 Myr before the SFR peaks). Dot- 
ted line shows the result if all the gas were consumed, /burst = /gas- Solid 
line shows the approximate model for the efficiency of disk survival in gas- 
rich mergers from| Hopkins et al!]j2009b} , roughly /burst ~ /gas (1 — /gas)- 
Small diamonds compare the effective EOS simulations used to analyze 
disk survival therein. Large circles compare the simulations here. Values 
/burst < /gas indicate disk survival and re-formation. Disk survival is com- 
parably efficient in full feedback and EOS models. 
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Figure 7. Optical image (as Fig.[2j of the stars in the relaxed merger rem- 
nant of the Sbc f merger (viewed edge-on and face-on). The violently re- 
laxed stars form a very extended halo, but this is mostly at surface densities 
~ 100 times smaller than the disk. The merger produces an SO/a (depend- 
ing on viewing angle), despite the fact that it just experienced an equal-mass 
(1:1) merger with no new gas accretion: almost any realistic situation will 
give even larger disks. 



rotational support following Hopkins et al. (2009b), by plotting the 
distribution of baryonic mass as a function of its angular momen- 
tum j/ j c , where j — |r x v| is the specific angular momentum and 
j c — r(GM to t(< r)/r) 1 / 2 is the angular momentum of a circular 
orbit for each particle. The distribution is strongly bimodal, with 
a peak near zero net angular momentum (bulge stars), and a peak 
near j/ j c = l, i.e. in a circular orbit. Fitting to this distribution (ei- 
ther cutting at the minimum between the two peaks or assuming the 
"bulge" peak is symmetric about j — and taking the "disk" mate- 
rial as the residual material: see Hopk ins et aLp 009b ) we can iden- 
tify the disk material as having an approximate rotational velocity 
of ^ 40kms _1 and mass of ~ 5 x 10 8 Mq, in good agreement with 
the observed baryonic Tully-Fisher relation (Bell & de Jong 2001 ; 
McGaugh 2005). This is not to say our simulations uniquely predict 
this relation - it follows from the initial disks lying on the observed 
relation by construction - but it does say that the disk is a "real" 
disk in as much as its progenitors were. Fig. [7] shows the optical 
image of the Sbc f merger relic, which also re-forms a prominent 
disk - based on the B/T ratio and presence of spiral arms, the face- 
on image is a typical Sa/b galaxy (edge-on SO/a). 

The key question is the quantitative efficiency of survival. 
Since the behavior of the dissipationless stars (i.e. how much they 
do or do not avoid violent relaxation) is identical in all models, we 
focus on the amount of gas which is able to survive the merger to re- 
form a disk. Hopkins et al. ( 2009b) parameterize this by examining 
the efficiency of gas consumption in the merger-induced starbursts. 
Gas which is not consumed (and not completely expelled by winds) 
will, inevitably, re-form a disk with whatever net angular momen- 
tum it carries (since it can cool). If all the gas were efficiently 
torqued (or locally forced to collapse in the merger), it would all 
contribute to the starburst, /burst = /gas - there would be no "resid- 
ual" disk gas survival. Hopkins et al. (2009b) argue, however, that 
at high gas fractions these torques become inefficient so the burst 
fraction scales sub-linearly with gas fraction as ~ / g as(l — /gas); 
this is the origin of the enhanced survival in gas-rich mergers. We 
compare the explicit feedback models to the EOS models used to 
calibrate this relation, and find they are statistically nearly identical. 
The starburst efficiency is clearly well below /burst = / gas ; i.e. disk 
survival is efficient. If anything there is slightly more disk survival 
in these runs, compared to the EOS runs, owing to the contribution 
of winds expelling some of the mass that would have contributed 
to the starburst bulge. 
© 0000 RAS, MNRAS 000, 000-000 



4 STAR FORMATION HISTORIES 

4.1 Starbursts 

Figure [8] shows the galaxy-integrated star formation histories of 
each simulation (averaged over « lOMyr intervals). The qualita- 
tive character of the models is similar, regardless of feedback de- 
tails. In prograde (e) orbits, resonant interactions lead to starbursts 
on first passage; however since the models here have some pre- 
existing bulge, these are not always strong (Barnes & Hernquist 
|1991|[T996] ); they are less prominent in the retrograde (f) case be- 
cause the systems are out-of-resonance. In all cases there is a star- 
burst on final coalescence, as the gas is channeled to the galaxy 
center. 

The integrated mass in stars formed in total and in the major 
starburst(s) are similar; as suggested by Fig. [5] However, the exact 
maximum SFR and shape of the SFH can differ significantly. It is 
already known that the duration, amplitude, and variability within 
the starburst are quite sensitive to the initial conditions, and treat- 
ment of sub-grid star formation recipe applied to the simulations 
(see |Cox et al.|2006c| , so this is not surprising. 

In all explicit feedback models, there is increasing variability 
in the SFR, owing to SF being concentrated in resolved GMCs and 
clusters, and feedback being associated with individual star clusters 
(and SNe) hence more stochastic in time and space. Interestingly, 
this suppresses some of the differences between prograde and retro- 
grade encounters in the weaker first-passage bursts; although quan- 
tifying this effect in detail requires a larger sample of orbits. 

In the low-mass (SMC and Sbc) cases, the starbursts are much 
more pronounced with explicit feedback models; this is because 
the EOS models, even with effective sound speed ~ lOkrns -1 , are 
sufficient to highly pressurize the systems and suppress a strong 
burst. In this case the allowance for molecular cooling makes the 
gas more compressible. 

However, especially in the higher-mass cases, expulsion of gas 
from the nuclei in feedback-driven winds (in the explicit feedback 
models) tends to "spread" the starburst in time. In particular, it gives 
rise to a long tail of SF as much of the material is not entirely 
unbound, but kicked into a fountain or stirred up within the disk 
and then re-condenses. The difference between this and the EOS 
case is especially clear in the prograde (e) models; there, with an 
EOS model, the gas is efficiently consumed in the starburst, leav- 
ing very little continuing SF afterwards. However, in the explicit 
feedback models, the SFR declines from a similar peak much more 
slowly, remaining at more than an order-of-magnitude higher SFR 
for ^Gyrs. 

4.2 The Kennicutt-Schmidt Law 

Figure [9] plots the positions of the simulations with explicit feed- 
back on the Schmidt-Kennicutt relation during the merger (the EOS 
models are placed on the relation by construction, and so are not in- 
teresting here). The agreement between these simulations and the 
observed relation for both isolated galaxies (the low-S gas regime, 
corresponding to when the galaxies are well- separated disks) and 
observed merging ULIRG/starburst systems is excellent. 

We remind the reader that there is no "tunable" parameter in 
these models - the feedback strength is not adjusted to match the 
relation, and the instantaneous star formation efficiency we assign 
in self-gravitating gas is unity, not the ^ 1 — 2% which is needed to 
match the Kennicutt relation. As we showed in Paper I and Paper 
II for isolated disks, the self-regulation of the SFR in these systems 
is a consequence of feedback; Paper I and Hopkins et al. ( 2 01 la) ) 
show that their location on the Kennicutt relation has almost noth- 
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Figure 8. Star formation histories in some representative merger simula- 
tions, in EOS and explicit feedback models. In all cases, prograde (e) mod- 
els produce more pronounced starbursts; the total mass consumed in the 
bursts is similar as well. The explicit feedback models predict more vari- 
able/bursty SFRs, because the ISM is less homogeneous. In low-mass sys- 
tems, the bursts are much stronger with explicit feedback models, because 
the EOS models heavily pressurize the gas relative to the disk V c (making 
it less compressible). At higher masses, the effect on the peak/duration of 
the starburst is more sensitive to other details. The "tail" of post-starburst 
SF is significantly enhanced in the explicit feedback models, especially in 
prograde mergers: in the EOS cases the starburst gas is all exhausted effi- 
ciently, but with winds present, some of this is expelled into a fountain that 
continuously recycles gas for the next couple Gyr. Starburst winds there- 
fore make merger-induced quenching/gas exhaustion (without AGN) less, 
not more, efficient. 



ing to do with the actual star formation prescription, once feedback 
is included (but is sensitive to the strength of feedback, if we were 
to artificially adjust it). Also note that there is no evidence for bi- 
modality in the relation, as suggested by Genzel et al. (2010); Daddi 
|et al.| ( [20T0] ); we discuss this in §[6] In Paper I, we show that the 
observed "cutoff" appears at sufficiently low E, but here it is less 
obvious since we always average over the half-SFR radius. 

4.3 Where Do Stars Form? 

In Figure[l0] we plot the radii that enclose 50% (R 5 o) and 90% (R 9 o) 
of the star formation in the merging galaxies, as a function of time[j 
At early stages the radii are large - basically the values for the iso- 
lated disks in seen in Paper II. At first passage, strong torques lead 
to shocks, dissipation, and rapid inflow; R50 drops rapidly as the gas 
piles up in a ~kpc- scale starburst. After this starburst, R50 increases 
again, until final coalescence when new inflows are driven to small 
scales and most of the star formation comes from a sub-kpc nuclear 
starburst. This is the standard scenario of merger-driven starbursts 
in simplified EOS models as well ( Hernquist 1989 ; Barn es & Hern-| 
|quist|1991||Mihos & Hernquist|1994b||1996| ). 

The evolution of R90 is significantly different. At early times, 
it also reflects the isolated disk sizes. However, going into first pas- 
sage (especially just before the first-passage starburst), and between 
this starburst and final coalescence, the radii actually become much 
larger, ^ 10 — 100 kpc. This is a consequence of star formation in 
tidal arms and the bridge region where the gas shocks between the 
two disks (especially prominent in prograde orbits). This is shown 
in Fig.[TT] where we plot the morphology of the young stars for the 
same simulations. In these stages, ~ 20 — 50% of the total SFR is 
contributed by this tidal/shock region. In contrast, it is well-known 
that in models with a strong effective EOS, large density contrasts 
in shocks are poorly captured and so the SFR in these regions tends 
to be small, - 10% of the total (see |Barnes 12004} |Karl et al.|2010) . 
Even though strong feedback is present, the gas in these simulations 
is (at least instantaneously) more compressible, so these shocks 
give rise to large density jumps that promote rapid star formation, 
significantly enhancing the contribution of the tidal region to the 
total SFR (seen in e.g. the no-feedback results with varied cooling 
prescriptions in Sait oh et al.|20 09). This may be further enhanced 
at even higher resolution, with the ability to better resolve strong 
density gradients in shocks on sub-pc scales; but at that scale, the 
effect of magnetic fields on shock compression is also important. 

5 MODELS WITH WEAK OR NO FEEDBACK 

In Paper I and Paper II, we discuss the properties of isolated simula- 
tions with no feedback. The models undergo catastrophic fragmen- 
tation with runaway gas collapse to arbitrarily high densities, turn- 
ing most of the gas into stars in a single dynamical time. We have 
not focused on these models here, because they are clearly ruled out 
by e.g. the observed Kennicutt (1998) relation. However, we briefly 
consider the differences between the predictions of the models here 
and a model with "weak" feedback implemented. Specifically, we 
consider an effective EOS model with q eos = 0, i.e. a cooling floor 
at 10 4 K or minimum turbulent velocity dispersion of lOkms . 
This is commonly adopted in historical models. This is not a no- 
feedback model; so it is not necessarily the case that the gas will 



6 When the galaxies are well-separated, we should consider the SFR within 
each separately; because our code is quasi-Lagrangian, we simply separate 
the particles belonging to each progenitor disk at all times, and take the 
appropriate radii separately (then plot the average of the two). 
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Figure 9. Kennicutt- Schmidt relation for the merger simulations with explicit feedback models. We compare the observed starburst galaxies in |Kennicutt| 
( |1998) and updated with high-redshift galaxies b y|Genzel et al.|{2010} (red x 's); shaded region shows the 90% range at each E gas from the compilations in 
those works and Bigiel et al. ( 2008); Daddi et al. ( 2010). Simulations are plotted at uniformly sampled times and with values measured at the half-SFR radius 
(averaged over 100 projections); before coalescence each merging galaxy is treated separately (populating the lower range). EOS models are constructed to 
lie on the relation by design, but the explicit-feedback merger models have a local star formation efficiency of unity in self-gravitating dense gas. Without 
feedback, the mergers lie a factor of ~ 20 — 100 above the observed relation (Paper I). With feedback - with no adjusted parameters - the KS-law is well- 
matched. As shown in Paper I and Paper II this emerges generically from feedback self -regulating the SFR (independent of the small-scale SF law); with the 
mechanisms included here, this extends even to extreme starbursts. There is no bimodality in the relation, though it may shift towards the upper envelope at 
high SFR. 



experience runaway fragmentation. In the isolated MW or SMC 
models, this is comparable to the observed velocity dispersion and 
that predicted from our full feedback model (see Paper II). How- 
ever, in more extreme systems - the mergers here, and even the 
isolated HiZ models - this is smaller than the predicted dispersions 
and the turbulent support needed to resist runaway gas fragmenta- 
tion. 

Fig.[l2]plots the SFR versus time (as Fig. [8), the half-SFR ra- 
dius versus time (as Fig. [ToJ, and the Kennicutt- Schmidt relation 
(as Fig. [9} for this model and our standard full feedback model. 
To highlight the differences, we focus on the most extreme case: 
the HiZ model and e orbit. This reaches the highest gas densities, 
SFRs, and gas densities of our simulations, so is most sensitive to 
the strength of feedback. With the weak feedback model, runaway 
collapse cannot be prevented. Early in the simulation (even before 
the galaxies are fully interacting, and especially on first passage), 
the gas undergoes catastrophic fragmentation and local clumps col- 
lapse to the highest resolvable densities (~ 10 6 cm -3 ), turning into 
stars rapidly. Because of this, the SFR is much higher on early pas- 
sages than final coalescence. But even in that coalescence, gas col- 
lapses rapidly as soon as it collects in filaments, before it can reach 
the nucleus. The effective radius of SF in this stage is ~ 5 — lOkpc, 
rather than the < 1 kpc traditional "nuclear starburst" we see in the 
simulations with full feedback; stellar clumps form during the in- 



flow and then sink collisionlessly into the bulgejjAs a result the 
SFR is more flat in the final coalescence (the burst enhancement 
is suppressed). This also means that the high nuclear gas densi- 
ties (Ega S ^> 10 3 Mqpc -2 ) observed in ULIRGs are not actually 
reached. The predicted SFR densities on the Kennicutt (1998} re- 
lation are uniformly an order of magnitude larger than observed 
(recall, the observations include high-redshift systems analogous 
to the models here). Since the initial effective EOS model is tuned 
to lie on this relation, this produces an apparent sharp bimodality 
in the relation. Finally, since there is no super-wind, the SFR de- 
clines rapidly post-merger (having exhausted all the gas and with 
little material at large radii to "rain back in"). At the end of the sim- 
ulation, we measure the efficiency of disk survival as in § |3.3| the 
surviving/re-formed disk mass in the weak-feedback simulation is 
only rsj 1 /3 of the mass of the re-formed disk in the full feedback 
simulation. 



7 Note that even during the final starburst, the total extent of the star- 
forming region in the full feedback models is quite large, > lOpc; quan- 
titative comparisons between similar simulations and observations of high- 
redshift extreme starbursts (sub-millimeter galaxies) shows that their spa- 
tial distributions agree well (see e.g. Younger et al. 2008b Narayanan et al. 
|2009||20T0||l"vison et al.|2009) . 
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Figure 10. Radial distribution of SF as a function of merger phase. For clar- 
ity we focus on one example (MW e), but the others are qualitatively similar. 
We plot the radii enclosing 50% (R 50 ) and 90% (R 90 ) of the total SFR@. In 
both EOS and full feedback models, the half-SFR radius initially reflects the 
isolated disks, but drops to ~kpc due to inflows induced by global torques 
on first passage, then to sub-kpc scales at the peak of the final coalescence 
starburst. The "compression" at the center is slightly larger in the explicit 
feedback models because of cooling to low temperatures allowing higher 
densities. The outer extent of SF (Rgo) is generally similar, except leading 
into first passage and between passages; here, more extended SF in tidal 
shocks (in arms and the inter-galaxy bridge region) is captured in the full 
feedback models, contributing ~ 20 - 40% of the total SFR. 
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Figure 11. Morphology of the extended SF plotted in Fig. [To] for the 
full feedback (left) and EOS (right) models. We plot only the very young 
(< lOMyr old) stars (projected surface density increasing with brightness 
with « 2dex stretch). We focus on the same example as Fig.[To|(MW e), at 
two times: just after first passage (top) and near final coalescence (bottom). 
The extended SF is much more prominent in the full feedback models, es- 
pecially after first passage in the tail/bridge regions. At final coalescence 
the kpc-scale starburst, containing most of the SF, has similar density and 
morphology in both runs, but is surrounded by more extended emission and 
a more prominent surviving star-forming disk in the full feedback model. 
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Figure 12. Comparison of a gas-rich merger with our full feedback model 
to one with artificially weak stellar feedback. We compare our standard HiZ 
e model to a simulation with identical initial conditions but an effective EOS 
with g e os = (cooling floor at 10 4 K). Top: SF history as Fig. [8] Without 
strong feedback, most of the gas is consumed by a massive burst within 
the disks when they begin to interact. The burst "enhancement" in final 
coalescence is much weaker. The post-merger SFR declines much more 
rapidly as the gas is exhausted, with no feedback-driven fountain to main- 
tain gas infall. The surviving disk fraction is lower by a factor ~ 3 in the 
weak-feedback model. Middle: Half-SFR radius as Fig.[l0| Even on coales- 
cence, most of the SFR in the weak-feedback case comes from large radii 
> 5 — 10 kpc, as inflowing gas undergoes catastrophic fragmentation and 
turns into stars before it reaches the nucleus. With full feedback present, the 
dip in R50 at coalescence reflects inflows driving a starburst (as Fig. [To) . 
Bottom: Kennicutt- Schmidt law as Fig. [9] With weak feedback, the highest 
observed surface densities (seen in merging nuclei) are never reached (gas 
fragments and turns to stars before reaching the central regions). The sepa- 
rated galaxies, constructed to lie on the KS law when isolated, also jump to 
much higher SFRs, predicting SFR densities an order of magnitude larger 
than observed and an apparently bimodal relation. 
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6 DISCUSSION 

We have studied major galaxy mergers with explicit models for 
stellar feedback that can follow the formation and destruction of 
individual GMCs and star clusters. Our models include star forma- 
tion only at extremely high densities inside GMCs, and the mass, 
momentum, and energy flux from SNe (Types I & II), stellar winds 
(both "fast" O-star winds and "slow" AGB winds), and radiation 
pressure (from UV/optical and multiply- scattered IR photons), as 
well as HII photoionization heating and molecular cooling. As a 
first study, we focus on simple, global properties, and compare 
them to those obtained from previous generations of simulations 
which did not follow these processes explicitly but instead adopted 
a simplified "effective equation of state" sub-grid model of the 
ISM. 

We find that most global consequences of mergers are similar, 
in models with explicit feedback and resolved ISM phase structure, 
and simplified EOS models. Some details, however, can be rather 
different. 

6.1 Galaxy Morphologies 

The relic structure and mass profile, and the presence of tidal tails 
in the merger, are predominantly determined by global gravitational 
torques, and so depend little on feedback details. This includes the 
integral properties (mass and effective radius) of the stellar relic of 
the kpc- scale starburst. 

Mih os & Hernquist| ( |1994a| ) argued that merger remnants 
are really two-component systems, and Ho pkins et al . (2008b a 
2009c a) later demonstrated this in comparison of observations (of 
ellipticals and recent merger remnants) and simulations with sub- 
grid feedback models. The outer profile is dominated by violently 
relaxed stars formed before the final galaxy merger, which are dis- 
sipationlessly scattered into the "envelope" at large radii. The inner 
portion is dominated by "starburst" stars which form via gas that is 
gravitationally torqued and dissipates, falls to the center, and turns 
into stars in a central starburst; this dominates the central ^kpc 
and extends inwards to form the central "cusps" in gas-rich merger 
remnants (see |Kormendy & Sanders |T992l |Hibbard & Yun|1999[ 
|Rothberg & Joseph||2004[ |Hopkins & Quataert||2011b| ); it is of- 
ten responsible for nuclear kinematically decoupled components 
Pernquist & Barnes|199"T||Hoffman et al.|2010| >. 

We find that this decomposition, and the behavior of the mass 
profile of each sub-component separately (as well as the total) is ro- 
bust to the inclusion of explicit or effective EOS models. This is re- 
assuring, since Hopkins et al. (2008b 2009a ) showed that the simu- 
lated mass distributions agree very well with observed mergers and 
spheroids, and |Cox et al.| ( |2006c| ) showed these results are robust to 
substantial changes in sub-grid parameterizations of feedback and 
star formation. The dissipationless profile should obviously be in- 
sensitive to the gas physics. The dissipative profiles are similar be- 
cause, as shown in Hopkins et al. (2008b ), their characteristic radii 
are set not by feedback "pressurizing" the material, but by the ra- 
dius at which the inflowing gas becomes self -gravitating, at which 
point it is no longer strongly torqued by the external stars and dark 
matter and can turn into stars relatively rapidly. The slope of the 
central mass profile is set by equilibrium between gravitationally 
driven inflow and star formation (Hopkins & Quataert 2011b). 

We may however expect differences in very low-mass, gas- 
rich mergers. In the EOS models, these are "pressurized" to the 
point where the effective sound speed is comparable to the circular 
velocity, and so tend to result in merger r emnants that are large r 
compared to true ellipticals at < 10 9 Mq ( [Hopkins et al.|2009a| . 



Allowing the gas to be more compressible (instantaneously), while 
retaining feedback, may resolve this discrepancy. 

6.2 Disk Survival 

Disk survival - the efficiency with which gas avoids consumption 
in the merger and can rapidly re-form a disk after - appears to be 
just as, if not slightly more, efficient in models with explicit feed- 
back (appropriately accounting for their properties at the time of 
merger). They follow the same scalings for the likelihood of disk 
survival (and gas angular momentum loss) as derived from simpli- 
fied models in |Hopkins e"taL](|2009b|). 

This is expected: |Hopkins et al.| ( |2009b| ) found disk survival 
did not depend on the sub-grid parameters of the EOS model (un- 
less feedback is very weak), and discuss the reasons for this. Gas 
is dissipative, so gas that survives the merger without turning into 
stars or losing most of its angular momentum will eventually re- 
form a rotating disk. The dominant torques driving angular momen- 
tum loss in mergers arise from resonant interactions exciting colli- 
sionless perturbations in the stars, which in turn force the gas into 
strong shocks ([Barne s & Hern quist|19^[l9^[Barnes|2002[[Hop^l 



|kins & Quataert|2010| 201 la| ). Angular momentum "cancellation,' 
viscosity, other hydrodynamic torques, and the direct torque from 
the presence of the secondary galaxy are typically much smaller 
effects. As a result, one can analytically derive a critical radius and 
resonant conditions that depend only on merger orbit, mass ratio, 
and galaxy structural properties (not the gas microphysics), outside 
of which gas will not be strongly torqued and will thus survive to 
re-form disks. With few stars, gas rich disks have weak collision- 
less torques, so give larger surviving disks. 

The same has also been seen in cosmological simulations with 
sub-grid feedback models ( |Governato et al.|[2009| |2010[ |Guedes| 
|et al.|2 011 ). Disks survive multiple major mergers by virtue of con- 
serving gas which immediately re-forms disks; this is enhanced by 
the presence of feedback-driven outflows which further suppress 
the mass turned into stars in the center. This is important for the ex- 
istence of realistic disks today ( Robertson & Bullock 2008 ; Hop- 
|kins et al.|[2009e[ |Stewart et al'||2009| >, and there appear to be a 
growing number of observed disks which have undergone recent 
merg ers ([Hammer et al.||2009b|a| |Yang et al.|[2009] |Puech et al.| 
[200^[Bundy et al.|2009||2010) . 

6.3 Starbursts and Star Formation 

The final starbursts are dominated by in situ star formation from gas 
which is torqued, falls into the nucleus and turns into stars at high 
densities. This is the "standard" scenario of merger-driven nuclear 
starbursts and star formation predicted in early simulations ([Barnes 
|& Hernquist|1991[[Mihos & Hernquist|1994b|[T996] ), and observed 
at the centers of nearby ULIRGs (strong, rapid inflows evident in 
their profiles of molecular gas, star formation, and metallicity; see 



e.g. |Tacconi et al.|[2002| |Sanders & Mir abel 1996, Kormen dy & 



Sanders|19921|Hibbard & Yun|1999[|Rothberg & Joseph|2004[|Ti- 
tus et al'|1997MReichardt et al.|20011|Michard 2006, F oster et al. 
Sanchez-Blaz quez et al.||2007| |Kewley et al.||2010[ [Rupke 
2010| |Soto & Martin|2010| ). |Cox et al.| < |2006c] showed that 
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et al. 



this is robust to large variations in the implementation of feedback 
in sub-grid models. As emphasized in many previous studies, we 
are not saying that most of the stars form in the starburst: this is typ- 
ically ~ 10% of the total stellar mass ( Jogee et al. 2009, Robaina 
|et al.|2009| [Hopkins et al.|2010[|Hopkins & Hernquist|2010| ), with 
the rest formed in the "isolated disk" mode (including times during 
the merger outside the "burst"). 
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Although gas clumps (into GMCs) as it flows in, feedback ef- 
ficiently disperses these GMCs after they turn just a few percent of 
their mass into stars. But the material is not expelled completely 
(which would shut down the burst entirely); most of the recycled 
gas is unbound locally from the GMC, but not globally from the 
galaxy. The typical velocities are the escape velocities from the 
GMC or local star cluster (~ lOkms -1 ) and so the gas is "stirred." 
This allows it to reach the nuclei in a coherent flow, so that gravita- 
tional torques will still tend to dominate. 

The more realistic feedback models, however, show signif- 
icantly higher star formation in tails and bridges (especially in 
shocks). This may resolve a long-standing discrepancy between 
merger models and observed star formation distributions in local 
mergers (Barnes 2004); a more quantitative comparison with mod- 
els intended to mimic these specific systems will be the subject of 
future work. 

In the explicit feedback models, the star formation is more 
strongly time variable (less "smoothed" by the effective EOS). This 
owes to the more inhomogeneous nature of the ISM - in extreme 
cases, a large fraction of the SFR can come from a few super star 
clusters. The burst "enhancement" (peak SFR) can be higher or 
lower, depending on the orbit and structural properties; as with 
sub-grid models, first passage bursts are much more sensitive to 
these details than final coalescence. In low-mass systems, it tends 
to be higher, because (as noted above) the EOS models may over- 
pressurize small disks by enforcing a minimum sound speed of 
~ lOkms -1 ; with molecular cooling the gas can more efficiently 
collapse to high densities in the nucleus. 

In all the explicit feedback cases, the "tails" of the starburst 
tend to be more broad; this owes to the effect of realistic stellar 
winds ejecting some material at modest velocities, so it falls back 
into the disk in a fountain and makes the starburst more extended 
(discussed below). 

6.4 The Kennicutt-Schmidt Relation 

The explicit feedback models predict SFRs that agree extremely 
well with observed merging galaxies on the Kennicutt-Schmidt re- 
lation. This is a major success of the models: unlike the EOS mod- 
els, where the kpc-scale SF law is imposed to match the Kennicutt 
law, the predicted Kennicutt relation in the explicit feedback cases 
is a true prediction (both in normalization, scatter, and power-law 
slope). Recall, star formation in these models occurs only in very 
dense, locally self-gravitating gas (scales <pc), with an instanta- 
neous local efficiency of unity - not the global ~ 1 % value of the 
Kennicutt relation. And without feedback (see Paper I & Paper II), 
the models predict a SFR far larger than observed, M* ~ Mg as /%n. 
With feedback included - with parameters not in any way adjusted 
to reproduce the observations but taken directly from stellar evolu- 
tion models, the observed relation naturally emerges. 

This is a consequence of feedback-regulated star formation: 
turbulent support in the gas dissipates locally on a crossing time, 
leading to collapse, which proceeds until a sufficient number of 
young stars form (regardless of how they form, in detail) to provide 
enough momentum flux in feedback to balance the dissipation. As 
such, we showed in Paper I that the SFR and Kennicutt relation are 
direct consequences of feedback, and are independent of the mi- 
crophysical star formation law once explicit feedback is included. 
We find the same here; we have re-run some of our simulations 
replacing our local self-gravity criterion for star formation with a 
simple density threshold with a very different efficiency (or chang- 
ing the small-scale density dependence), as discussed in § [2] and 
the predicted Kennicutt law is nearly identical. 



The same feedback mechanisms regulating star formation in 
isolated galaxies extend to the regime of extreme merger-induced 
starbursts. This is no mean feat: the typical surface densities and 
star formation rates are orders of magnitude larger than the isolated 
disk counterparts of these galaxies. The gas surface densities in the 
bursts here extend to ~ 100 — 1000 times larger than the mean sur- 
face density of local GMCs ! We are therefore well into the regime 
where the medium is entirely molecular and high-density (although 
it may still exhibit phase structure). In Paper II, we showed that 
SNe and other "thermal" feedback mechanisms become unimpor- 
tant in gas above densities ~ 1 — 10 cm -3 (typical of the diffuse 
gas in GMCs), because the cooling times are extremely short com- 
pared to the dynamical times. We therefore expect that radiation 
pressure is the most important feedback mechanism in these star- 
burst regions, as we found there for the most dense regions within 
"normal" GMCs. It is a remarkable success of the models, though, 
that they interpolate naturally into even this extreme regime, with- 
out requiring any adjustment. In contrast, most analytic models for 
the Kennicutt-Schmidt relation have been forced to assume some 
ad hoc break or transition in the scaling at high densities. 

We see no strong bimodality in the predicted relation. But 
this is not necessarily inconsistent with recent observational claims 
along these lines ( [Genzel et"aT|2 010 ; Daddi et al. 2010); our merger 
and isolated galaxy simulations overlap closely with the obser- 
vations from those papers. What we see is a continuum between 
these systems, which is completely consistent with the observa- 
tions (see e.g. Narayanan et al. 2011a). The apparent bimodality 
stems from a strict choice of two separate CO conversion factors; 
in future work, we will incorporate detailed models for molecu- 
lar emission to make a quantitative comparison with these observa- 
tions. Narayanan et al. ( 201 lb) find that when applying a functional 
form for the CO-H2 conversion factor which depends on the phys- 
ical conditions in the ISM to the observed galaxies in Figure [9] a 
slope of rsj 1.8 emerges. Our results are quite consistent with those, 
and we find a slope of ~ 1.7; this suggests that stars do form some- 
what more "efficiently" (in terms of consumption time) in systems 
like mergers which occupy the high- surface density region of the 
Kennicutt-Schmidt plot, though mergers do not lie on a different 
track of the relation than disks. The increased efficiency likely owes 
to a larger fraction of gas in the dense phase, a result which has been 
noted in observations of local ULIRGs (e.g. | Juneau et al.|2 009). 

6.5 Comparison to Models with Weak Feedback 

Recently, there have been some different conclusions reached in 
the literature from the study of high-resolution simulations with 
very different implementations of stellar feedback. |Teyssier et al.| 
( 2010) argued that merger-induced starbursts might not occur "in 
situ" from inflows, but instead from catastrophic fragmentation of 
tidal arms, bars, and other instabilities. The fragments would turn 
a large fraction of their mass into stars very rapidly, and then sink 
as supermassive stellar clusters to the galaxy center via dynamical 
friction, building the bulge. They also suggest a sharp bimodality 
in the Kennicutt-Schmidt relation, with the mergers and isolated 
disks being well-separated, as a consequence of this induced frag- 
mentation in the mergers (though they do discuss how this might 
be smoothed out). Bour naud et al.| p010b| ) argued that disk sur- 
vival is very inefficient in their simulations, for similar reasons - 
the gas undergoes dramatic fragmentation and "clumps" scatter off 
one another, providing new torques. 

We believe that all of these differences stem from the fact that 
these simulations include very weak stellar feedback, which is un- 
able to disperse dense gas clouds (GMCs and especially clouds in 
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the dense tidal or starburst regions of mergers) after they form. In 
fact, we reproduce all of the above results if we simply re-run our 
mergers with an effective EOS model with much weaker "feed- 
back strength," q eos = 0. Like the q eos = models, the Teyssier 
|et al.| ( [20T0] ) models adopt an effective EOS, but one in which the 
"median" ISM temperature is much colder than the sub-grid sound 
speeds used in the EOS models here, which makes catastrophic 
fragmentation possible. Without strong feedback, the clumps can- 
not then disperse and recycle their mass. Specifically, the models 
either do not include explicit feedback mechanisms or include only 
SNe; however, we show in Paper I and Paper II that SNe have al- 
most no effect on the dispersal of dense gas, because the cooling 
time in that gas is > 100 times shorter than the dynamical time. 
As a result, the SNe can stir up the "diffuse" medium but cannot 
recycle dense gas. In our simulations, the most dense clouds are 
destroyed by a combination of radiation pressure in both UV and 
IR, photo-ionization, and momentum from O-star winds. 

In models that allow cooling but do not include these explicit 
feedback mechanisms, therefore, gas cools and collapses into dense 
GMC-like objects, which steadily contract and form stars, but are 
not efficiently "re-mixed" into the ISM on a short timescale. Col- 
lapse (hence star formation) in the dense regions runs away, giving 
rise to the dramatically enhanced "fragmentation mode" of star for- 
mation and apparent bimodality in the Kennicutt relation. In Paper 
II we show that even if slow star formation is forced in the dense 
regions, the lack of mixing allows the dense regions to move bal- 
listically, with relatively small interaction cross-sections and long 
lifetimes. In this limit, gas is effectively no longer collisional! All 
the dissipation occurs within individual clumps, while the relative 
motions between clumps cannot be dissipated, just like relative mo- 
tions of stars; the gas behaves like a collection of "super-massive" 
star or dark matter particles, and loses orbital angular momentum 
only via processes much slower than shocks, such as scattering and 
dynamical friction. 

In contrast observed GMCs are short-lived, with lifetimes of 
a few free-fall times (few Myr), and turn just a few percent of 
their mass into stars befor e dispersing ([Z uckerman & Evans|1974{ 
IWilliams & McKee|[T997l|Evans|19991|Evans et al.|2009b) . This 
appears to hold even for most of the most massive (10 — 10 9 Mq) 
GMC complexes observed in local mergers (Rand & Kulkarni 
[T990l|Plane"sas et al.| 1991] | Wilson et al.|2003||2006) . So even if 
most of the gas at a given instant is in dense sub-units (GMCs), the 
gas will be recycled through the diffuse ISM on a short timescale, 
and therefore experiences normal hydrodynamic forces. So long as 
the GMC lifetime is short compared to the total duration of the 
strong torques in the merger (~ 500 Myr), this will be true. 

Likewise, since inflows in the final coalescence of a major 
merger require only a couple of dynamical times to reach the cen- 
ter, for fragmentation to "beat" inflow as a starburst driver the frag- 
menting clumps in the inflow would have to have a very high in- 
stantaneous star formation efficiency, M* ~ M gas /%n. This hap- 
pens without feedback; however, the observed Kennicutt relation 
in ULIRGs and other mergers (even high-redshift systems) implies 
an efficiency a factor of ^ 50 smaller. We note that direct observa- 
tions of ULIRGs find that the final coalescence star formation rate 
is indeed dominated by a nuclear starburst (see references above). 

A critical conclusion is that, in any model which resolves the 
formation of dense clouds, it is equally important to include and 
resolve the physics that may destroy those clouds. 
© 0000 RAS, MNRAS 000, 000-000 



6.6 Future Directions 

The results here are a preliminary study of how realistic stellar 
feedback and ISM structure influence the results of merger sim- 
ulations. We have found that the most basic global integral prop- 
erties of star formation and galaxy morphology tend to be similar 
to those inferred in previous studies with a simplified treatment of 
the gas physics. However, we also find that more detailed structural 
and kinematic properties are sensitive to the gas physics. We do not 
expect this to be the case outside of the central ^kpc, where the 
relic galaxy is dominated by stars formed before the merger and 
violently relaxed. But kinematic substructures in the central region 
may be sensitive to how the gas collapses and turns into stars. It will 
also be particularly interesting to examine the effects of more real- 
istic star formation and enrichment models that can separately fol- 
low stellar winds and SNe Types I & II on the age, metallicity, and 
abundance (a/Fe) gradients in the galaxies; these can place strong 
constraints on the merger history and role of dissipation in galaxy 
formation (see e.g. Torrey et al .|2012| and references therein). 

In a companion paper, we examine the star clusters formed in 
these simulations. The mass/luminosity distribution, spatial loca- 
tions, formation time distribution, and physical properties of these 
clusters represent a powerful constraint on small-scale models of 
the ISM and star formation. 

We have also restricted our focus to major mergers. Studies 
of mergers with varying mass ratios suggest that the qualitative be- 
haviors discussed here should not depend on mass ratio for ratios 
to about 3:1 or 4:1, and even at lower mass ratios they can be con- 
sidered similar but with an "efficiency" of inducing starbursts and 
violent relaxation that scales approximately linearly with mass ra- 
tio pernquist & Mihos|| 19951 |Cox et"aT][2008l [Naab & Burkert| 
|2003[ | Younger et a~ 2008a). But at small mass ratios the dominant 
role of mergers may be more subtle: inducing resonant disk per- 
turbations ( [D'Onghia et al.|2009||20T0|) and disk heating ([Purcell 
et al.|20091|Stewart et al. |2008|| Walker et al.|1996l|Hopkins etaL 
2008c i Moster et al. 2010 ), for which the gas response may depend 
significantly on its phase structure. 

We note that recent studies comparing cosmological sim- 
ulations done with GADGET and the new moving mesh code 
AREPO ( Springel 2010) have called into question the reliability 
of smoothed particle hydrodynamics for some problems related to 
galaxy formation in a cosmological context ( Vogelsberger et al. 
1 201 l||Sijacki et al.|201 l||Keres et al.|201 l[|Bauer & Springel|201 1[ 
Torr ey et al.||2011| ). However, we have also performed idealized 
simulations of mergers between individual galaxies and found ex- 
cellent agreement between GADGET and AREPO for e.g. gas-inflow 
rates, star formation histories, and the mass in the ensuing star- 
bursts, for reasons that will be discussed in (Hayward et al. 2 01 lb| ). 
The discrepancies above are also minimized when the flows of in- 
terest are supersonic (as opposed to sub-sonic), which is very much 
the case here ( |Kitsionas et al.|2009||Frice & Federrath|2010[|Bauer| 
|& Springel|20TT) . 

Our new models allow us to follow the structure of the gas 
in the central regions of starburst systems at high resolution. This 
makes them an ideal laboratory for studying feedback physics un- 
der extreme conditions in, say, the center of Arp 220 and other very 
dense galaxies. In another companion paper, we examine the prop- 
erties of the large starburst winds driven by stellar feedback in the 
mergers, which can reach ^ 10 — 500MQyr -1 outflow rates and 
so have potentially major implications for metal enrichment and 
self-regulation of galaxy growth. We have also here, for clarity, ne- 
glected AGN feedback in these models, but we expect it may have 
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a significant effect on the systems and their outflows after the final 
coalescence. For example, it may strongly suppress the otherwise 
quite high post-merger SFRs we see in these simulations, which - 
without something to rapidly "quench" them - may make it difficult 
or impossible to explain the observed abundance of high-redshift, 
low-SFR galaxies. With high-resolution models that include the 
phase structure of the ISM, it becomes meaningful to include much 
more explicit physical models for AGN feedback. 
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APPENDIX A: SIMULATION IMAGES 

Figures |A1|A8| present the images of the rest of our simula- 
tions, at different viewing angles and times during the merg- 
ers, in the style of Figures |1|2| Movies and additional im- 
ages are available online (jhttps : / /www, c fa. harvard. | 
|edu/~phopkins/Site/Movies_sbw_mgr . html). 
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Figure Al. Morphology of the e (prograde) merger of the Sbc galaxy (a dwarf starburst). We show two perpendicular viewing angles (left and right panel 
sets); each shows the optical (as Fig. [2] left) and gas (as Fig.[T] right). The time since the start of the simulation (in Gyr) is labeled. From top to bottom, the 
merger stages are (1) pre-merger infall, (2) just after first passage, (3) apocenter/turnaround, (4) final merger, (5) just after nuclear coalescence. 
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Figure A3. As Fig.|AT] for the SMC (SMC -mass dwarf) e merger. 
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Figure A4. As Fig. |Al| for the SMC f merger. 
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Figure A5. As Fig. |Al| for the MW e merger. 

© 0000 RAS, MNRAS 000, 000-000 



Hopkins et al. 



t = 0.5 

m 

10 kpc 
t = 0.7 


30 kpc 


t = 0.5 

10 kpc 
t = 0.7 


30 kpc 


10 kpc 


J* 

30 kpc 


10 kpc 


30 kpc 


t= 1.3 

m 

10 kpc 


r \ 

30 kpc 


t= 1.3 
10 kpc 


30 kpc 


t = 1.5 

lb 

10 kpc 


30 kpc 


t= 1.5 
10 kpc 


30 kpc 


t = 1.8 

* 

10 kpc 


V 

30 kpc 


t= 1.8 

Hp 

10 kpc 


30 kpc 


t = 2.5 

9 

10 kpc 
t = 3.3 


"'V'"- 

* > 

30 kpc 


t = 2.5 

1 

10 kpc 
t = 3.3 


* 

30 kpc 


10 kpc 


i 

30 kpc 


10 kpc 


30 kpc 



Figure A6. As Fig.|AT] for the MW f merger. 
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